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We study parametric amplification of Kaluza-Klein (KK) modes in a higher D-dimensional 
generalized Kaluza-Klein theory, which was originally considered by Mukohyama in the narrow 
resonance case. It was suggested that KK modes can be enhanced by an oscillation of a scale of 
compactification by the d-dimensional sphere S d (d = D — 4) and by the direct product 5"* 1 x 
S 2 (di + d2 = D — 4). We extend this past work to the more general case where initial values of 
the scale of compactification and the quantum number of the angular momentum I of KK modes 
are not small. We perform analytic approaches based on the Mathieu equation as well as numerical 
calculations, and find that the expansion of the universe rapidly makes the KK field deviate from 
instability bands. As a result, KK modes are not enhanced sufficiently in an expanding universe in 
' these two classes of models. 
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I. INTRODUCTION 



Recently, much interest has been focused on higher-dimensional theories. The original idea goes back to Kaluza 

ON . 

, and Klein 111] , who considered that the unification of the interactions in Nature may be realized in the gravitational 

o . 

theory in five dimensions. At present, there are many higher dimensional theories which are the generalizations of the 
Kaluza-Klein theory, for example, D = 10 superstring theory and D = 11 M-theory j^j. Very recently, the universe 



of brane models W has received much attention as a solution to the hierarchy problem between the weak and Planck 



Oh, scale. In order to hide extra dimensions, we generally assume that they are within a very small compact internal 
<D . 

space. Since the typical scale of the internal space is the Planck length l p \ ~ 10 _33 cm, which is much smaller than the 

> ■ 

size of the external space, we can reduce higher-dimensional theories to the familiar four-dimensional gravitational 

& '• 

5r , theories. 

There are various ways for Kaluza-Klein reductions to four dimensions. As was pointed out in Ref . |^] , it is possible 
to judge whether such reductions are appropriate from a cosmological point of view. For example, in the model where 
the geometry is described by a direct product of the four-dimensional Minkowski spacetime M A and 1-sphere S 1 
considered by Kolb and Slansky 0], particles of Kaluza-Klein (KK) modes which are so called pyrgon are produced. 
Since the energy density of this massive particle decreases as pkk ~ a ~ 3 5 this density will dominate over the energy 
density of the universe in the radiation dominant era if pyrgons are overproduced in the early stage of the universe 



*electronic address:shinji@gravity.phys. waseda.ac.jp 



1 



and do not decay by entropy productions. Mukohyama |5[ extended this model, and considered the compactifications 
both by a d-dimensional sphere S d with d — D — 4 and by a direct product S dl x S d2 with c?i + di = D — 4. In 
these classes of models, since the reduced effective potential in the four-dimensional theory lacks a global minimum, 
we need to introduce the Casimir effect as a one-loop quantum correction to give a stable ground state. Then a 
scalar field er which corresponds to the scale of the compactification oscillates around the minimum of its potential, 
and is finally trapped at the ground state a = 0. Since KK modes acquire masses which are expressed by the 
oscillating er field, we can expect the amplification of KK modes by parametric resonance. This picture is similar to 
the scenario of preheating after inflation where scalar fields coupled to an inflaton field are resonantly enhanced in 
the oscillating stage of inflaton |^,[| . In the case of a chaotic inflation model with a quadratic potential, particles are 
most efficiently produced in the broad resonance regimes where the amplitude of inflaton <f> is initially large and the 
coupling between cf> and a coupled scalar field x is strong ||. Since the expansion of the universe gradually reduces 
the amplitude of inflaton, this makes the x field jump over many instability and stability bands. The x field can 
grow quasiexponentially, overcoming the diluting effect by the cosmic expansion. Finally the field reaches so called 
narrow resonance regimes, after which particle creation terminates. In the case that the x field is initially located 
in narrow resonance regimes, it is known that particle creation is inefficient because the field soon deviates from 
instability bands by the cosmic expansion. As for the excitation of KK modes in the models of compactifications by 
S d and S dl x S d2 , since the a field behaves as a massive scalar field around the minimum of its potential, we can 
apply analytic approaches based on the Mathieu equation fllpf which have been much investigated in the context of 
preheating at the linear stage of the system ||[ll) . However, since the investigation in Ref. § about the enhancement 
of KK modes is limited in narrow resonance regimes in which parametric resonance is weak, we will extend this work 
in the more general range of parameters. It is of interest whether catastrophic particle production will occur or not in 
the case where the value of a is initially large and the quantum number of the angular momentum I of KK modes is not 
small. If KK modes are strongly enhanced by parametric resonance, this would result in some important cosmological 
implications. They will become good candidates for dark matter in the case that these excited particles are suitably 
diluted by entropy productions. 

This paper is organized as follows. In the next section, we perform compactifications by a <i-dimensional sphere S d 
and by a direct product S dl x S d2 . The Casimir energy is introduced to stabilize the scale of the internal space. In 
Sec. Ill, we study the background evolution of the scale factor and the a field, and examine the structure of resonance 
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for KK modes. We will investigate whether the enhancement of KK modes efficiently occurs or not both by analytic 
approaches and numerical integrations. We present our conclusion and discussion in the final section. 

II. THE MODEL 



We consider a model in D = d + 4 dimensions with a cosmological constant A and a single scalar field 



S = I d D x^7j 



; R~2h--g MN d_ 



M<PON<, 



(2.1) 



,2k 2 " 2 C 

where g~MN and R 2 /8ir = G are the Z?-dimensional metric and gravitational constant, respectively. R is a scalar 
curvature with respect to §mn- The third term of the action ( |2.l| ) denotes a Klein-Gordon action in D dimensions. 
We first compactify extra dimensions to a d-dimensional sphere S d . Then the metric §mn can be expressed as 



ds 2 D = gMN<lx M dx N — g yLV dx^dx v + b 2 ds 2 , 



(2.2) 



where g^ is a four-dimensional metric, 6 is a scale of the d-dimensional sphere whose present value is b$, and ds 2 d is 
a line element of the d-unit sphere. We expand the scalar field </> on the sphere as 



= b n d/2 ^ <f>imY l 

I .in 



(d) 

1.711 ' 



(2.3) 



where Y$ is the spherical harmonics on the d-sphere with positive integer I and a set of d — 1 numbers m. Then we 



obtain the reduced action in four dimensions, 

d 



S= I d 4 



1 

2k 



,,R + d(d 



l,m 



l(l + d-l) , 2 
2b 2 ^ 



(2.4) 



where k 2 /8tt = S 2 /(87rVj) is Newton's gravitational constant with V® the present volume of the internal space. R is 
a scalar curvature related with g^ v . Since this does not take the ordinary form of the Einstein- Hilbcrt action because 
of a time dependent factor (b/b ) d , we perform a conformal transformation as follows 



g^ = exp 

where a is the so-called dilaton field which is defined by 

a = oqitl 

CO 



d— I g^, 



(2.5) 



~d(d + 2) 
2k 2 



1/2 



(2.6) 
(2.7) 



Then the four-dimensional action in the Einstein frame can be described as 



S = / (Txy^g 



2k 2 2" ^ " iV ' 2 

i 



(2. 



where R is a scalar curvature with respect to . U\{a) is a potential of the a field which we will give a specific form 
below, and the mass Mi{a) of the <pi m field is expressed as 

M f{a) = liJ ^"^ e-^. (2.9) 

As was mentioned in Ref. Q, the potential Ui(a) does not have a local minimum if we do not introduce the Casimir 
effect as a one-loop quantum correction. Taking into account this effect, U\{a) can be written in the form 



1 _ e -2(d+2)(r/ CT o i e -d<r/(To _ d + 4 c -(d+2)cr/cr 



d+2 d+2 



(2.10) 



J7i(a) = a 
with 

a = JdTW - (2 - n) 

The first, second, and third terms in ( |2.1C| ) appear due to the Casimir energy, the cosmological constant, and the 
curvature of the internal space, respectively. The potential Ui(a) acquires a local minimum at a — for the case of 



d > 2 because of the presence of the first term in (|2.10| ). We depict Ui(a) in Fig. 1 for the case of d = 2. It has a local 



maximum at cr*(> 0), which depends on the extra dimension d. a*/ao decreases with the increase of d. For example, 



c*/(To = 0.50 for d — 2, and <t*/<7o = 0.22 for d = 6. Since Co becomes larger with the increase of d [See Eq. (2.7)], 
<7* itself does not necessarily decrease with the increase of d. For example, <r* = 0.20rn p i for d = 2, and er* = 0.21m p i 
for d = 6. In order to result in the final state a — which corresponds to the present value b = 6 0j initial values of 
a are required to be < <7j < c* (we assume <7j > 0), where the subscript i denotes the initial value. Then a evolves 
toward the minimum of its potential, and begins to oscillate around a = 0. In this oscillating stage, we are concerned 
with whether KK modes are effectively enhanced or not by parametric resonance. 

In the case that the compactification is done on a direct product S dl x S d2 of the e?i-dimensional sphere S dl and 
^-dimensional sphere S d ' 2 with d\ + d 2 = D — 4, the final four-dimensional action and the potential U\{a) take 
different forms from the S d case. The procedure is the same as in the S d case. First, we consider the D-dimensional 
metric 

ds 2 D = g M Ndx M dx N = g^dx^dx" + b\ds 2 di + b\ds\, (2.12) 



where b\ and b 2 are scales of the sphere S dl and S d2 , ds 2 di and ds 2 d2 are line elements of the di-unit sphere and c^-unit 
sphere, respectively. 



Expanding the <j) field in the D-dimensional action (2.1 



as 



,-d 1 /2,-d 2 /2 
"10 u 20 



V 61 1 Y {dl) Y {d2) 

l 1 ,l 2 ,m 1 ,m 2 



(2.13) 



where 610 and 620 are present values of 61 and b 2 , VJ^I and are the spherical harmonics on the di-sphere and 

d2-sphere, respectively, we obtain the four-dimensional action similar to the action (2.4) in the S d case (see Ref. || 
for details). Next, we perform a conformal transformation 

' d\(j\ d 2 a 2 



g^ v = exp 

where Oiii = 1,2) and Uio are defined by 



<?10 020 



en = a l0 hi 



0iO 



ho J ' 

+ 2)1 1/2 



2k 2 



Taking into account the Casimir effect as in the S d case and imposing the condition 



01 _ 02 _ 

010 020 00 



with 



cr = 



(di +d 2 )(di + d 2 + 2) 



2k 2 



the final expression of the four-dimensional action yields 



S = I d A Xy/—g 



9 X! {g^d^hhmnnzdv&hhmn^ + Mf 1 l 2 {^)4 l 'i 1 l 2 m 1 , 



Ii,l2,mi,m2 

where the potential of the a field is 



Ui(a) = a 



di+d 2 + 2 



e -2(d 1 +d 2 +2)a/a„ , e -(d 1 +d 2 )a /a _ d l + d 2 + ^ c -(d 1 +d 2 +2)a / a„ 

d x + d 2 + 2 



with 



a 



d!+d 2 + 2 
2{di +d 2 + 4)k 2 



di(di-l) , d 2 (d 2 -l) 



bio 



bio 



(2.14) 



(2.15) 
(2.16) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



The mass of Mi 1 i 2 (a) of the 4>i 1 i 2 m 1 m 2 ne ld is given by 



liih + di-V, l 2 {l 2 + d 2 -\) 
h 2 h 2 

"10 u 20 



e -{d 1 +d 2 +2)a-/a-Q ^ 22) 



The terms in the square bracket of Eq. (2.2C) correspond to changing d in that of Eq. ( p. 1C| ) to d\ + d 2 . Hence the 
shape of the potential U\(a) in the S dl x S d2 case is the same as in the S d case. U\(a) has a local minimum at a = 
when either of d\ and d 2 is greater than 1 . In the next section, the dynamics of a and </>/ m fields are studied in two 
classes of models of (2.8) and ( 2.19| ). 



III. THE DYNAMICS OF THE SYSTEM 

In this section, we investigate the excitation of KK modes due to the oscillating scalar field a by parametric 
resonance. We assume that the four-dimensional spacetime and the a field are spatially homogeneous, and adopt the 
flat Friedmann-Robertson- Walker metric as the background spacetime 

ds\ = g^dx^dx" = -dt 2 + a 2 (i)dx 2 , (3.1) 

where t is the cosmic time and a(t) is the scale factor. The background equations for a and a can be expressed by 
the action (fj|) or (|2.19| ) as 









K 2 " 




3 




3^-f 








a 



2 <J Z + U 1 (a) 



(3.2) 



a) = 0, (3.3) 

where we have neglected the contributions from KK modes. However, in the case that KK modes are enhanced 
significantly, this affects the evolution of background equations. We briefly comment on this back reaction issue at 
the end of this section. 

As we have already mentioned, the initial value of a is required to be smaller than c* in order not to lead to an 

expansion of the internal space. Let us consider the case where a is initially located close to cr, with <Ji < a*. In this 

case, defining a dimensionless parameter 

_f(d + 2)^ (S d case), 

£ ~ i (di + d 2 + 2)^ (S dl x S d2 case), ^ A > 

it usually exceeds of the order of unity at the beginning with d > 2 or d\ + d 2 > 2. For example, in the S d case, 
e(U) = 2.0 when d = 2, and e(^) = 1.8 when d = 6. This is the situation which was not considered in Ref. ||. After 



U{'(0) 



-y/2(d-l) 
bo 



di(di-l) , d 2 {d 2 -l) 



(3.5) 



a short stage of rolling down, the a field begins to oscillate around the minimum of its potential at a = 0. After that, 
the a field behaves as a massive scalar field whose mass is given by 

(S d case) 

rf7+3J 77-^ ~ 7^ x ^ 2 case ) 

The evolution of this field is the same as a massive inflaton field (f> with a chaotic potential V((^>) = m^0 2 /2 in the 
reheating phase. Making use of the time-averaged relation 

1 



we easily find from Eqs. J3.2| ) and ( |3.3| ) that the evolution of a and a can be approximately written as 

2/3 



(3.6) 



(3.7) 



a(t) /<jq = a(t) cosra CT i, 



(3.8) 



where t co denotes the time when the coherent oscillation of a starts, and a(t) is the dimensionless amplitude of the 



a field which decreases as a(t) ~ 1/t by the cosmic expansion. Eq. (3.7) shows that the universe evolves as matter 
dominant when the a field oscillates coherently. Hence the energy density of a decreases as p a ~ 1/t 3 . 

Let us study the excitation of KK modes during the oscillating stage of the a field. For each compactification of 
the S d and S dl x S d2 cases, we expand 4>i m {S d case) and <fii 1 i 2 m 1 m 2 (S dl x S d2 case) fields by the Fourier component 



as 



4>hhm im2 = ( 2?r )3/2 /"( a fe0W2mim 2 (*) e + a l^W 2TOim2 (*)e* k ' X ) 



(3.9) 



(3.10) 



where at and at are the annihilation and creation operators, respectively. Then we find that the temporary parts 
(j)f m (t) and 0£/. 2mi77l2 (*) obey the following equations of motion: 

~k 2 l(l + d-l) 



A k 4. -l^j,* 



b 2 



-(d+2)a/a 



bf m = (S d case), 



(3.11) 



hjh+dx - 1) l 2 {h + d 2 -l) \ {dl+d2+2 )*/* 

"10 U 20 



= (S' dl x S"* 2 case). 



(3.12) 



Hereafter, we express both fields <fif m and ^ 1 i 2 m 1 m 2 as exce Pt the case that distinctions of both fields are required. 
Defining a new scalar field tp k = a 3 / 2 (j) k , Eqs. ( 3.1l| ) and (3.12) are written as 



0fc +w k <p k = 0, 



(3.13) 



where for the S d case, 



1(1 + d- 1) (d+2W<70 _ 3 /2o ff 
6§ 



+ ■■ 
4 V a a 



(3.14) 



and for the x S"* 2 case, 



k 2 (h(h+di-l) l 2 (h+d 2 -l)] _r^_L^^„/^ 3 /2a d 2 



6 2 



b 2 20 



p -(d 1 +d 2 +2)a/a a — 

4 \ a a 2 



(3.15) 



As for the initial conditions of the ip k field, we choose the state that corresponds to the conformal vacuum as 

1 



<Pk(ti 



-, Vktti) = -iu k (ti)ip k (ti). 



(3.16) 



y/2u k (ti 

The second terms in Eqs. ( ^.14| ) and ( 3.15} ) are resonance terms which would lead to the amplification of KK modes. 
When e defined by Eq. (3.4) is greater than of the order of unity, resonance terms are suppressed and do not play 
relevant roles for the enhancement of the ip k field. However, e decreases under unity and a begins to oscillate coherently, 



we can expect the parametric resonance phenomenon. In this case, Eq. (3.13) is reduced to the Mathieu equation, 

d 2 

— 1 ip k + (A k -2qcos2z)ip k = 0, (3.17) 

where for the S d case, 

21(1 + d- I) A (k/m a ) 2 



A, = 



d- 1 



l(l + d- 1). , „._. . 
g= d _ 1 ' (d + 2)a(t) 7 



(3.18) 



(3.19) 



and for the S dl x S d2 case, 



. 2(d 1 + d 2 ){h(h + di - 1) + h(h + d 2 - l)(b w /b 20 ) 2 } (k/m a ) 2 

A k = 7T-T-, — T-T-, -CT. ~ h 4- 



q = 



d 1 (d 1 -l) + d 2 (d2-l)(ho/b2 ) 2 
eh + d 2 ){h(h +d 1 -l) + l 2 (l 2 + d 2 - l)(b 10 /b 20 ) 2 } 

d 1 (d 1 -i) + d 2 (d 2 -i)(b w /b 20 y 



(di + d 2 + 2)&(t). 



(3.20) 



(3.21) 



z is defined as z = rat. We have neglected the contribution of the third term in Eq. ( 3.14 ) and ( 3.15 ) since this term 
is not important in the oscillating stage of the a field M . 



We can analyze the excitation of KK modes by making use of the stability and instability chart of the Mathieu 



equation |10|. Parametric resonance occurs when KK modes stay in the instability regions sketched in Fig. 2. These 
regions are typically described by narrow resonance regimes of q < 1 and broad resonance regimes of q 3> 1. As is 
found in Fig. 2, the instability bands are broader for larger values of q as long as A k is not so large. Generally, particle 
creation in narrow resonance regimes is not so efficient as in the case of broad resonance regimes [9 . Since the past 
work about the excitation of KK modes S is limited in narrow resonance regimes q < 1, we extend this work in the 
more general case with q 1. However, note that the results obtained in the context of preheating with a quadratic 
potential || are not directly applied to the present case, because the relation of A k and q is different. 
First, let us first examine the S d case. As is found in Eq. ( |3.18 ), Ak approaches the constant value 

_ 21(1 + d-l) 

B= (d-1) ' (3 - 22) 

with the passage of time for fixed values of d and I. Since B is greater than unity for the positive integer d and I 
with d > 2, parametric resonance does not occur in the first resonance band around Ak — 1 as was pointed out in 
Ref. ||. For example, when d — 2 and I = 1, B = 4. The initial value of q (= q co ) at which the analysis based on the 
Mathieu equation becomes valid depends on the value of a when the a field begins to oscillate coherently. Numerical 
calculations show that this corresponds to e(t co ) sa 0.4, after which the evolution of the a field can be described as 
Eq. ( |3l| ). This means that the (d + 2)d(t) term in Eq. ( |3. 19| ) and the (d 1 +d 2 + 2)a(t) term in Eq. ( |3.21| ) are restricted 
to be smaller than unity. In order to realize the situation g> 1, we need to choose larger values of I, or l\ and l 2 . 

Let us investigate two cases of q co < 1 and q co S> 1. First, we consider the case of q co < 1 with the extra dimension 
d = 2. When I = 1, since A k = 4 + 4(fc/m (T a) 2 and q co = 8a(t co ) m 0.8 by Eqs. (|3.18|) and J3.19|) , the tp k field is 
initially located in the second instability band around Ak = 4 for low momentum modes. However, the value of q 
becomes smaller with the decrease of the amplitude a(t) as a(t) ~ \ jt by the cosmic expansion. In Fig. 3, we depict 
the evolution of the a field in this case. We choose a(ti)/ao = 0.4 as the initial value of the a field. After the first 
stage of rolling down, a begins to oscillate with the initial amplitude &(t co ) ~ 0.1. The coherent oscillation starts by 
the same amplitude as long as the initial value of a is located in the region of 0.1 < u(ti)/ao < ov/co ~ 0.5. After 
several oscillations, the amplitude drops down a(t) < 0.02, which corresponds to q < 0.1. Although the iff. field stays 
in the instability band since A k approaches A k — > 4, the growth rate fi k of the ip k field becomes smaller with the 
decrease of q. In the j-th resonance band with A k = j 2 , q <C 2j 3 / 2 , and j > 2, fi k is expressed as pi 
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where 8 changes in the interval [0, ir/2}. In the second instability band with q < q co m 0.8, fik cx q 2 . The growth rate 
(ik is initially of the order 0.01, and after that [ik rapidly decreases. The growth rate of KK modes <pk — ar z / 2 yk is 
approximately expressed as 

d(pk _-3/2 ( ^ H 



where we have used the relation tp k ~ e MfcZ . Since the Hubble expansion rate at t = i co is -ff/m ~ 0.28, which exceeds 
the growth rate of tpk, 4>k decreases at the beginning. Moreover, since fik and H decrease as \ik cx q 2 oc 1/t 2 and 
H cx 1/t, the increase of <pk can not be expected in the second instability band of the narrow resonance regime q < 1. 
We have numerically confirmed that <fik does not increase in the case of d = 2 and Z = 1 (see Fig. 4). If there is a set 
of d and I which satisfies B = j 2 with positive integer j > 3, the narrow instability bands exist around the regions of 
Ak = j 2 - However, parametric resonance is much more inefficient than in the second resonance case, since the growth 



rate jiu decreases with the increase of j for q < 1 as is found by Eq. (3.23). Even if d and / are changed, the first 



resonance does not occur because B defined in Eq. (3.22) is always greater than unity. In the narrow resonance of 
q < 1 with the S d compactification, we can conclude that KK modes are not enhanced sufficiently in any values of 
the extra dimension d and the positive integer I. 

Consider the case of q ^> 1 with the S d compactification. As we have already mentioned, larger values of / are 
required to lead to the condition g> 1. In order to judge whether parametric resonance is efficient or not, we should 



know the relation between Ak and q. In the S case, we obtain the following relation from Eqs. ( 3.18 ) and ( 3.19 ) 



2 (fe/m CT ) 2 

A ^idrm) q+4 -^- (3 - 25) 



The efficiency of resonance strongly depends on the value of (d + 2)a(t) in Eq. (3.25). We have numerically confirmed 
that this value is approximately (d + 2)a(t co ) « 0.4 when the a field begins to oscillate coherently for any values 
of d. This corresponds the relation of Ak ~ 5g co + 4(fc/m cr a) 2 at the beginning of the oscillating stage. Although 
there are many instability bands in the regions of q 3> 1 which are generally called broad resonance bands, these 
bands become rare in the regions of Ak > 5q (see Fig. 2). Moreover, the tangent of the line A k = 2q/{(d+ 2)a(t)} 
gets larger with the passage of time because of the decrease of the amplitude a(t) as a(t) ~ 1/t, which means that 
parametric resonance becomes inefficient. Even in the case that the <fk field is initially in an instability band, the 
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expansion of the universe makes the field deviate from the instability band. In the model of a massive inflaton with 
a coupled scalar field x m preheating after inflation, the \ ne ld moves on the path of Ak = 2q + k 2 / (to 2 a 2 ) with the 
decrease of q as q ~ 1/i 2 . In this case, there are many instability bands as well as stability bands in the regions of 
Ak > 2q, the x field stochastically increases when passing through instability bands. Since the tangent of the A^-q 
path does not change, x particles can be efficiently produced overcoming the diluting effect by the cosmic expansion 
P). In the present model, the excitation of KK modes is very weak because the regions of Ak > 2q/{(d + 2)er(i)} are 
mostly occupied by stability bands. Although the tpk field passes through instability bands, these are very few in the 
regions of Ak > 2q/{(d + 2)a(t)}, and parametric resonance is not efficient. Even the low momentum modes close to 
the line Ak — 2q/{(d+ 2)a(t)} are not enhanced sufficiently As time passes, both Ak and q decrease by the cosmic 
expansion, and finally approaches Ak = B = 21(1 + d — l)/(d — 1) and q = 0. In Fig. 5, we show the evolution of the 
KK mode tpk with d = 2, I = 100, and k = 0. In this case, q co as 8000. We find that the KK mode does not grow as 
in the case of q < 1. Namely, the growth rate fj,k of the ifk field does not surpass the Hubble expansion rate even in 
the case of q 3> 1. As we have already mentioned, the initial relation of Ak and q hardly depends on the values of 
o~(t CQ ) and d, which means that parametric resonance is also inefficient even if we choose larger values of cr(U) and d. 
Next, let us consider the compactification by S dl x S d2 . In this case, since the value 



_ 2(di + d 2 ){l 1 (l l +d 1 -l) + l 2 (l 2 + d 2 - 1)( W& 20 ) 2 } 

d 1 (d 1 -i) + d 2 (d 2 -i)(b 10 /b 20 y 



B = TU TTT^m TT7T. 71 \2 ( 3 ' 26 ) 



is greater than unity for positive integers l\, l 2 and d\, d 2 with d± > d 2 and d\ + d 2 = D — 4, the first resonance does 
not occur as in the case of S d [^2|. When B is written by a positive integer j(> 2) as B — j 2 and the value of t^o 
is smaller than unity, the ifk field is located in an instability band at the initial stage. However, since q decreases 
as q ~ 1/t by the cosmic expansion and the growth rate of KK modes is very small for j > 2, the excitation of KK 
modes is inefficient in narrow resonance regimes q < 1. In the case of q ^> 1, since the relation of Ak and q is written 
by Eqs. ( fT2C| ) and ( ^2l| ) as 



Ak H — i j , n \~(+\ q + 4 2 ' (3.27) 

(di + d 2 + 2)a(t) a 1 

the discussions in the S d case can be applied by changing d\ + d 2 to d. We have numerically confirmed that the a 
field begins to oscillate coherently at (di +d 2 + 2)a(t co ) w 0.4 for any values of a(U), and d\, d 2 . Moreover, since a(t) 
decreases as a(t) ~ 1/t, parametric resonance is also ineffective as in the S d case. 

We conclude that catastrophic particle creation does not occur in two classes of models by S d and S dl x S d2 
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compactifications in any values of parameters. As a result, the back reaction effect by KK modes is not important in 
these models. 

IV. CONCLUDING REMARKS AND DISCUSSIONS 

In this paper, we have studied the excitement of Kaluza-Klein (KK) modes in a higher Z?-dimensional generalized 
Kaluza-Klein theory. We have considered two classes of models where the extra dimensions are compactified on the 
sphere S d with d — D — 4 and the direct product S dl x S d2 with d\ + d 2 = D — 4. Such compactifications give rise to 
a term which originates from the curvature of the internal space in the potential U\(a) of a dilaton field a. Since this 
potential does not have a local minimum to stabilize the scale of the internal space, we introduce the Casimir effect 
due to a one-loop quantum correction. Then the a field oscillates as a massive scalar field around the local minimum 
at a = 0, which corresponds to the present value of the scale of compactifications. 

Since the KK field acquires a mass which can be expressed by the a field by compactifications on S and S dl x S d2 , 
we can expect that KK modes will be enhanced by parametric resonance. The past work on this issue [^) is restricted 
in the case of narrow resonance regimes where the resonance parameter q is smaller than unity. However, in the 
similar situation of preheating after inflation with a quadratic potential, it is well known that resonance with q 1 
is much more efficient than in the case of q < 1 . Hence we extend past work on the excitation of KK modes to the 
case of q S> 1 by making use of the stability and instability chart of the Mathieu equation. 

In the case of q < 1, the first resonance does not occur for both compactifications by S d and S dl x S d2 . Although 
there are some situations where the KK field is located in the second, third, • • • instability bands at the beginning of 
the coherent oscillation of the a field, parametric resonance soon becomes ineffective with the decrease of q due to the 
expansion of the universe. In this case, since the creation rate of KK modes can not surpass the Hubble expansion 
rate, the enhancement of KK modes is inefficient. 

Even in the case of g > 1, we have found that the growth of KK modes does not take place both by analytic 
approaches and numerical integrations. The a field begins to oscillate coherently when the amplitude of the a field 



drops down to e « 0.4, where e is defined by (3.4). Since the relation of resonance parameters A k and q are expressed 
as A k = 2q/{(d + 2)a(t)} + A(k/m a a) 2 for the S d case, and A k = 2q/{(d 1 +d 2 + 2)d{t)} + A(k/m a a) 2 for the S dl x S d2 
case, the KK field exists in the region of A k > 5q where instability bands are few at the beginning. The amplitude 
of a decreases with the passage of time, and the KK field evolves in the regions where instability bands are further 
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few. As a result, KK modes are not relevantly enhanced even for the case of q ^> 1. We have numerically confirmed 
this fact, and found that the excitation of KK modes is inefficient in any parameters in two classes of models of 
compactifications. 

Since wc find that KK modes are not overproduced by parametric resonance by compactifications of S d and S dl x S d2 , 
this kind of compactification may not be ruled out from a cosmological point of view, because the energy density 
of KK modes will not overdose the universe in the radiation dominant era. However, at the stage of preheating 
after inflation, scalar fields coupled to an inflaton field can be strongly enhanced by parametric resonance |1 . If the 
KK field is coupled to inflaton, the enhancement of KK modes would also occur in the preheating stage. Recently, 
Mazumdar and Mendes jl3| considered the excitement of the dilaton field a as well as the Brans-Dicke field during the 
preheating phase in generalized Einstein theories. Since they compactified the extra dimensions on torus which does 
not have the curvature of the internal space, the potential of dilaton U\(a) does not appear in their model. Taking 
into account the growth of metric perturbations during preheating pi], it was found that the dilaton field a can be 
effectively enhanced even when dilaton does not couple to inflaton. Although they did not consider the enhancement 
of KK modes, there will be a possibility that KK modes are strongly enhanced in the preheating phase even in the 
case where the KK field does not directly couple to inflaton by the growth of metric perturbations. In this model, 
back reaction effects would play an important role for the termination of resonance. Although it is technically difficult 
to deal with back reaction issues including second order metric perturbations in a consistent way, it is of interest how 
KK modes are enhanced in the preheating phase. These issues are under consideration. 
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Figure Captions 

FIG. 1: The potential U\(o) which is obtained by introducing the Casimir effect in 

the S d compactification with d = 2. The potential has a minimum at a = 
and a local maximum at er*/(7o = 0.50. If we choose larger values of d, <T*/ao 
becomes smaller. The shape of this potential in the case of the S dl x S d2 
compactification is the same as in the case of the S d compactification. 

FIG. 2: The schematic diagram of the Mathieu chart. The lined regions denote the 

instability bands. There exists narrow instability bands around A k = j 2 and 
q < 1 with positive integer j. Although there are many instability bands in the 
regions of q 3> 1, they are few for Ak > 5q. 

FIG. 3: The evolution of a as a function of t in the case of the S d compactification 

with d = 2, I = 1 and fc = 0. We choose the initial value of a as a/<r = 0.4. 
After the first stage of rolling down, the a field begins to oscillate coherently as 
Eq. (3.8). The dimcnsionlcss amplitude a(t) in Eq. (3.8) decreases as a(t) <~ 1/t 
with the initial value of a{t co ) w 0.1. 

FIG. 4: The evolution of the real part of the Kaluza-Klein mode (f>k as a function of t 

in the case of the S d compactification with d = 2, I = 1, and k = 0. Although 
the (f>k field is initially in the second instability band with q < 1, the expansion 
of the universe makes parametric resonance ineffective. As a result, the growth 
of 4>k can not be expected. 

FIG. 5: The evolution of the real part of the Kaluza-Klein mode 4>k as a function of t in 

the case of the S d compactification with d = 2, I = 100, and k — 0. Although 
the initial value of q is large as q ^> 1, the <j>k field moves in the regions of 
Ak > 2q/{(d + 2)a(t)} where instability bands are few. Hence parametric 
resonance is inefficient. 
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